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EFFECTS OF AIRFOIL PROFILE ON THE TWO-DIMENSIONAL 
FLUTTER DERIVATIVES FOR WINGS OSCILLATING 
IN PITCH AT HIGH SUBSONIC SPEEDS 


By John A. Wyss and James C. Monfort 


SUMMARY 


Aerodynamic lift and moment flutter derivatives were determined &t 
high subsonic speeds for a series of two-dimensional airfoils varying in 
thickness and thickness distribution. The wings were sinusoidally oscil- 
lated about the quarter-chord axis at Mach numbers from about 0.5 to 0.9. 
The corresponding reduced frequency ranges varied from 0.045 to 0.45 
at M = 0.5 and from 0.025 to 0.25 at M=0.9. An evaluation of the 
results indicated that wing profile and angle of attack have major 
effects on the flutter derivatives at speeds exceeding the Mach number 
for steady-state lift divergence. In general, at supercritical Mach 
numbers the trends of the magnitudes of the oscillatory lift coefficients 
were qualitatively indicated by the trends of the nonoscillatory coeffi- 
cients, with phase angles, except for the l2-percent-thick airfoil, 
having only moderate deviation from subsonic theory. The variations in 
the magnitude of the moment derivative and in its phese angle, resulted 
in a trend toward instability at supercritical Mach numbers. In particu- 
lar, for airfoils of equal thickness the effect of an extreme forward 
location of maximum thickness was destabilizing in that negative -aero- 
dynamic damping existed, implying the possibility of a single degree of 
freedom type of flutter. Decreasing airfoil thickness delayed the large 
deviation from subsonic theory to higher Mach numbers. 


INTRODUCTION 


This report is concerned with the evaluation of the effects of air- 
foil profile on the lift and moment flutter derivatives as measured, by 
means of pressure cells, on harmonically vibrating two-dimensional wings 
at high subsonic speeds. It is well-known that theory does not account 
properly for such factors as flow separation and shock formation, hence, 
the aircraft designer must of necessity look to experimental values 
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whenever such mixed-flow conditions are encountered. Numerous previous 
investigations at lower speeds, such as those by Clevenson and Widmayer 
(ref. 1) and by Halfmen (ref. 2), may be cited. With the use of a dif- 
ferent measuring technique, the present work extends these previous 
investigations to higher Mach numbers so that emphasis may be placed 


upon supereritical speeds for which information is meager or nonexistent, 


Since wing profile may be expected to have а significant effect on 
mixed-flow conditions, several models were used to determine the effects 
of wing thickness and thickness distribution on the flutter derivatives. 
NACA 65۵ series symmetrical airfoils, 12, 8, and 4 percent thick, were 
used along with two other 8-percent- thick airfoils with their maximun 


thickness at about 16 апа 63 percent of the wing chord. The models were | 


oscillated about the quarter-chorá axis at Mach numbers from 0.5 to 0.9 
with reduced frequency ranges from 0.045 to 0.45 and from 0.025 to 0.25, 
respectively. Reynolds numbers, based on the airfoil chord, varied 
from 5 to 8 million. 


SYMBOLS 

a velocity of sound in undisturbed air, ft/sec 

b: wing semichord, ft 

сі. dynamic section lift coefficient 

Cm ` dynamic section moment coefficient about quarter point of 
chord 

f frequency of oscillstion, cps 

k reduced frequency, ae 

Y 

M Mach number, — 

Ма: -oscillatory aerodynamic section moment on wing about axis of 
rotation, positive with leading edge up 

Ра. oscillatory aerodynamic section lift on wing, خی ون‎ EM 
upwards 


free-stream dynamic pressure, 1b/sq ft 


ү. free-stream velocity, ft/sec 


b uli 21% 
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а, oscillatory angular displacement (pitch) about axis of rota~ 
tion, positive with leading edge up, radians 


Om mean angle of attack about which oscillation takes place, deg 
e phase angle between oscillatory moment and position с, posi- 
tive for moment leading a, deg 
0 phase angle between oscillatory lift and position a, posi- 
tive for lift leading a, deg 
ш circular frequency, 2xf, radians/sec 
d Реті 
Ln magnitude of dynamic lift-curve slope, M ,per radian 
da 2baa 
den е-19 
一 -一 magnitude of dynamic moment-curve slope, 一 三 一 一 | ,Per radian 
да, lib“ ас. 
а . : 
e sin 9 aerodynamic damping component in phase with angular velocity 
ей 


APPARATUS AND METHOD 


Models and Instrumentation 


The 12- and 8-percent-thick airfoils, NACA 65A012, 65A008, 2-008, 
and 8774008} profiles, were of wood-rib and wood-stressed-skin construc- 
tion built around steel spars at the quarter chord, which was the axis 
of rotation. Several wood spars at other chordwise locations were used 
to minimize spanwise twisting since the models were driven from one side. 
The h-percent-thick model, of МАСА 65A004 profile, was machined from 
solid aluminum with a parting line in the chord plane. The upper and 
lower halves of this model were bolted and doweled together. Each model 
had a chord of 2l inches and a span of 18-1 + inches. The gaps between 
the ends of the models and tunnel walls were sealed with sliding spring- 
loaded felt pads or brass strips which moved with the models. 


lAn NACA 8۱۷۴۵330 airfoil was modified to a symmetrical section by 
using the lower surface coordinates for both upper and lower surfaces 
and then reducing the thickness ratio to 8 percent. 
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In figure 1, the model profiles are illustrated to show the varia- 
tion of thickness and thickness distribution. The reference model, 
NACA 65A008, is marked to indicate the locations of the pressure cells. 
Model instrumentation consisted of 15 flush-type pressure cells (see 
refs. 3 and 1) and 15 pressure orifices along the midspan of each sur- 
face of each model. The pressure orifices adjacent to each pressure 
cell were used for two purposes: (1) as a means to determine the time- 
average chordwise pressure distribution with the use of a multiple 
mercury manometer, and (2) to provide an internal reference pressure 
for the pressure cells. The tubes from each cell and from the adjacent 
pressure orifice were interconnected at the manometer. In order that 
the internal reference pressure of the pressure cells would be essen- 
tially steady, about 50 feet of 1/16-inch tubing was used from the 
orifice to the manometer and back to the pressure cell. 


Two l4-channel oscillographs were used to record the instantaneous 
electrical difference of the output of each pair of cells (proportional 
to the pressure difference between the upper and lower surface at each 
chord station) and to record the summation of all cells (proportional 
to the variation of the lift force). The output of an NACA slide-wire 
position transducer, proportional to the model angle of attack, was 
simultaneously recorded. 


Tunnel, Model Drive System, and Tests 


. The models were oscillated in the two-dimensional test section in 
the Ames 16-foot high-speed wind tunnel (ref. 5). The two-dimensional 
channel was about 20 feet long and 16 feet high. A view of a model in 
place and a diagrammatic sketch of the drive system are presented in 
figure 2. The drive rods and sector arm attached to the model were con- 
tained within one of the channel walls. 


Records were obtained with Mach number and mean angle of attack 
constant for frequencies from № to 10 cycles per second at intervals of 
4 cycles per second and for an amplitude of +10. Date are presented 
for mean angles of attack of 0° and 2° and for Mach numbers from 0.5 to 
about 0.9. Sample oscillograph records which illustrate the necessity 
for harmonic analysis at the higher Mach numbers are given in figure 3. 
The lift was evaluated by a l2-point harmonic analysis of three consec- . 
utive cycles of the sum trace. The pitching moment was evaluated by a 
12-point harmonic analysis of the individual cell traces for one cycle. 


, Since the investigation was conducted in а closed-throat tunnel, 
the effects of wind-tunnel resonance must be accounted for either by 
avolding conditions in which tunnel-wall effects are significant or by 
correcting the results for the effects of the tunnel wells (refs. 6 ` 
and 7). Calculations made at the Langley and Ames Laboratories employing 
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the single-doublet-line, single-control-point solution described іп ref- 
erence 7 ylelded the following results for a tunnel height of 16 feet, 
wing chord of 2 feet, and Mach number of 0.7: At frequencies of 10, 20, 
and 40.66 cycles per second, the magnitudes of the coefficients were 
increased by 3.8, 5.0, and 4.7 percent, respectively, due to the presence 
of the tunnel walls. These results indicate that, for the conditions of 
the calculations, the effect of the tunnel walls was small. However, for 
mixed-flow conditions, the application of such corrections based on 
potential flow would be questionable; hence, to minimize tunnel-wall 
effects, all data obtained at frequencies within 10 percent of the tun- 
nel resonant frequency (refs. 6 and 7) have been omitted. Although the 
use of such a procedure does not mean tunnel-wall effects have been com- 
pletely eliminated over the entire frequency range, it is felt that 
tunnel-wall effects are not a predominant factor in the trends of the 
data. 


For a discussion of other factors influencing the precision of the 
data, the reader is referred to references 3 and L. 


RESULTS AND DISCUSSION 


A tabulation of the measured derivatives is contained in tables I, 
ІІ, III, IV, and V for the МАСА 65A012, 65А008, 65A004, 2-008, and 
877A008 airfoils, respectively. The results concerning lift derivatives 
are first discussed and are presented in figures 4 to 10, followed by a 
discussion and the presentation of the moment derivatives in figures 11 
to 15. КЫ | 


Lift 


Experimental values for the reference model for three representa- 
tive Mach numbers are presented in figure h as a function of reduced 
frequency. Tm this figure, as in subsequent figures, the absolute mag- 
nitude of the flutter derivative is expressed in terms of the slope of 
the lift curve per radian and the corresponding phase-angle relationship 
between the lift vector and model angle of attack in degrees. Theoreti- 
cal values at Mach numbers of 0.5, 0.6, and 0.7 may be obtained from the 
work of Dietze (refs. 8 and 9), and at Mach numbers of 0.8 and 1.0 from 
Minhinnick (ref. 10) and Nelson and Berman (ref. 11), respectively. 


In this figure it may be noted that at 0.49 and 0.79 Mach numbers 
the flutter derivatives tend to increase with increasing reduced fre- 
quency; furthermore, there seems to be a large variation in the phase 
angle at low values of reduced frequency at 0.79 Mach number. However, 
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Mach number appears to have had a greater effect on the data than did 
frequency at 0.91 Mach number. 


Typical results as a function of Mach number are presented in fig- 
ure 5 for the reference model, the NACA 65A008 airfoil. The lines show- 
ing the theoretical values are identified at one end by the frequency in 
cycles per second to which they pertain. Since theoretical values have 
been computed in the cited references only at certain Mach numbers which 
have already been indicated, an interpolation was necessary to obtain 
values at intermediate Mach numbers. Although such an interpolation | 
inherently involves some error, a consistent set of values was neverthe- 
less established and was used for the purpose of determining the effects 
of varying airfoil shape. 


To distinguish between the various frequencies, the experimental 
and theoretical values are each faired with the same type of line. For 
example, the experimental and theoretical values for a frequency of 8 
cycles per second are each shown with a solid line. Examination of the 
experimental data for & frequency of B cycles per second indicates that 
the trends of both experiment and theory were the same at low Mach num- 
bers. As Mach number increased, a large decrease in the magnitude of 
the experimental derivative occurred, accompanied by a variation of 
phase angle such that the trend toward increasing lag was reversed. 
Although the agreement with theory was not precise at the lower Mach num- 
bers, it may be seen that the general trends for all frequencies were 
nearly the same. 


The data from figure 5 are presented in a different form in fig- 
ure 6; the experimental magnitude has been divided by the theoretical ` 
magnitude, and the theoretical phase angle has been subtracted from the 
experimental phase angle. These quantities are also shown as a function 
of Mach number. If the experimental and theoretical values exactly 
agreed, the ratio of the magnitudes of the derivatives would be 1, while 
the difference in phase angle would be O. The faired lines represent 
the average deviation from theory for the entire frequency range up to ` 
40 cycles per second. f | m i 1 


It is of interest to note that the individual points do not indicate 
an entirely random scatter about the mean line for the various frequen- 
cles. For example, examination of the points for HO cycles per second 
in the top portion of the figure shows that these points are usually the 
uppermost value at each Mach number. Hence, this figure not only provides 
some indication of the range of the experimental values, but illustrates 
the fact that, although the values depend on frequency, the general varia- 
tions with Mach number are represented by the faired average curves. 


The use of the average deviation from theory appears to be justified 
since it is representative of each model. For example, in figure 6 it 
may be noted that all the experimental points lie within a comparatively 
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narrow band along the faired curves with the exception of the higher 
frequeneies in the upper portion of the figure. In fact, a band of 
width £0.15 in the upper portion of the figure апа a band of width +100 
in the lower portion of the figure would contain about 80 percent of all 
the experimental points. These results are typical of all the models. 
It might be noted that the averaging process used has the effect of 
removing frequency as a parameter. It should be noted that each model 
was oscillated at the same amplitude and through the same range of fre- 
queneies, hence the average deviation from iheory indicates the over-all 
effects of airfoil shape and the general trends of the data. 


Effect of thickness distribution.- The effects of the variation of 
thickness distribution as indicated by the curves showing the average 
deviation from theory over the frequency range tested are summarized in 
figure 7 for mean angles of attack of 09 and 2°. It would appear from 
this figure that the main effect of the chordwise location of maximum 
thickness was on the magnitudes of the derivatives rather than on phase 
angles, although no systematic trend is apparent. 


Effect of wing thickness.- The results showing the effects of wing 
thickness are presented in figure 8. At an angle of attack of 0°, wing 
thickness appears to have had a much more pronounced effect than wing- 
thickness distribution (fig. 8(а) as compared to fig. 7(а)). As might 
be expected, the primary effect of reducing wing thickness was to delay 
any large deviation from theory to a higher Mach number. 


At an angle of attack of 2° (fig. 8(b)), large differences over the 
entire range of Mach numbers occurred between the models in the magni- 
tudes of the derivatives. 


Comparison with steady-state results.- In order to examine whether 


any relation existed between unsteady and steady-state results, a com- 
parison with steady-state results obtained from the time-average chord- 
wise pressure distributions for mean angles of attack of O° and 2° is 
made in figures 9 and 10. In these figures, the steady-state data have 
been normalized with the Prandtl-Glauert value of the theoretical lift- 
curve slope. It may be recalled that the Prandtl-Glauert curve is also 
obtained as an end condition as the frequency of oscillation approaches 
Zero. 


Examination of these figures indicates that although there appears 
to be some parallelism or similarity between the steady and unsteady 
curves, the comparison between the steady and unsteady values is at best 
only qualitative. For example, in neither figure 9 nor figure 10 do the 
unsteady and steady-state curves coincide throughout the entire range of 
Mach numbers. It should also be noted that, with the exception of the 
NACA 65A012 airfoil at a mean angle of attack of 2° (fig. 10(b)), the 
unsteady values approached theory more closely than did the steady-state 
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values, particularly at the lower Mach numbers, that is, from M = 0.5 
to 0.7. Although the effect of the higher frequencies in increasing the 
level of the curves for the unsteady case may in part account for the 
differences between the curves, this effect is small. However, the one 
characteristic that is common to both the unsteady and steady curves in 
almost every case is a trend toward a reduction in magnitude at the 
highest Mach numbers. The Mach number at which this trend initiates 
cannot be precisely delimited, nevertheless, for the three NACA 65A-series 
airfoils at a mean angle of attack of 0° (fig. 10(а)), the unsteady lift 
trend appears to be associated with the steady-state flow changes which 
occur above the Mach number for lift divergence. 


It would therefore appear that as a first approximation the Mach 
number for lift divergence may be taken as a criterion for the onset of 
significant changes in the trends of the unsteady values, and that this 
trend toward a decrease in the magnitude of the unsteady values is 
related to the trend of the steady-state data. It should be pointed out 
that this conclusion is not as evident for the NACA 2-008 and 877A008 
airfoils (fig. 9) and for the NACA 65A004 airfoil at a mean angle of 
attack of 2° (fig. 10(b)), since these figures indicate that the correla- 
tion between the Mach number for lift divergence and the initiation of a 
downward trend of the unsteady values is not precise and they may differ 
by as much ав 0.1. However, it is felt that there is sufficient evidence е 
presented in figures 9 and 10 to indicate that steady-state values may 
prove useful as a qualitative indication of the trends of the unsteady- 
state coefficients at supercritical Mach numbers. 


For the steady-state condition the phase angle is, of course, zero; 
therefore no corollary for the phase angle with relation to the oscilla- 
tory condition is possible. However, except for the 12-percent-thick 
wing, the phase angle shows only a moderate deviation from theory through- 
out the speed range of the present investigation. 


Moment 


The moment derivatives for the reference model as a function of 
reduced frequency for several Mach numbers are presented in figure 11 
and as a function of Mach number in figure 12. A comparison of these 
figures indicates that even though there may have been a greater effect 
due to frequency on the moment derivatives than had been the case for 
the lift derivatives, from figure 12 it appears that the effects of Mach 
number are similar for all frequencies. Hence, the effects of airfoil 
profile are again compared on the basis of the faired average curves in 
figure 12 which represent the average deviation from theory over the 
entire frequency range. 
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In contrast to the lift results previously presented in figure 6, 
the magnitudes of the moment derivatives greatly exceeded the theoretical 
values, along with a much larger variation of phase angle as compared 
with theory. These results may be attributed to the fact that the com- 
parison is between very small quantities in regard to the magnitude of 
the derivatives, since the moment is taken about the quarter-chord axis, 
and to small movements of the center of pressure which would be reflected 
in large changes of phase angle. The general trends of the results, 
nevertheless, are represented by the faired average curves. 


Effect of thickness distribution.- The effects of the variation of 
the chordwise location of maximum thickness are shown in figure 13. An 
apparent characteristic of the NACA 2-008 airfoil, with a forward loca- 
tion of maximum thickness, isa large shift toward a lagging phase angle 
as Mach number increased above 0.8, such that the phase angle lagged 
theory by 80° and 90° at angles of attack of O° and 29, respectively. 
The effects of such large shifts in phase angle are discussed in rela- 
tion to subsequent figures. 


Effect of wing thickness.~ The effects of wing thickness on the 
moment derivatives are shown in figure 14. As might be expected, the 
primary effect of decreasing wing thickness was again to delay any large 
variations to a higher Mach number. 


Instability.~ Since there was such a large variation at the higher 
Mach numbers from the subsonic theoretical values, it is of basic impor- 
tance to examine the damping-moment derivatives directly to determine 
whether instability, or the existence of negative aerodynamic damping 
(implying the possibility of a single degree of freedom type of flutter), 
which is not predicted by the theory, existed at these speeds. The 
average damping-moment derivatives for the entire frequency range are 
therefore presented in figure 15. Also included in this figure are 
dashed lines indicating average values derived from theory for the corre- 
sponding frequency range. 


The effect of wing=-thickness distribution on aerodynamic damping is 
shown in figure 15(а) for each mean angle of attack. It may be noted 
that there was a trend toward instability for each model, with the 
NACA 2-008 airfoil becoming abruptly unstable at about 0.85 Mach number 
at 09 and 2° angles of attack. It would appear that stability about 
the quarter-chord axis increased as maximum thickness was moved toward 
the trailing edge. 


The effect of wing thickness on the aerodynamic damping moment is 
shown in figure 15(Ъ) for each angle of attack. Although the trend 
toward instability does not appear at 09 angle of attack for the 
NACA 65A004 profile, the susceptibility of the thinner wing to negative 
aerodynamic damping is clearly indicated at the 2° mean angle of attack. 


м 
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CONCLUSIONS 


Within the limitations of speed range and angle-of-attack variation 
of the investigation, the following general conclusions may be drawn: 


1. Section profile has a major effect on the flutter derivatives 
at speeds exceeding the Mach number for steady-state lift divergence. 


2. It appears that the variation in angle of attack has an effect. 
as important as the effect of the variation in profile. 


3. Іп general, at supercritical Mach numbers, s qualitative evalua- 
tion of the results indicated that the trends of the magnitudes of the 
oscillatory lift coefficients were indicated by the trends of the non- 
oscillatory lift coefficients, with phase angles, except for the 
12-percent-thick model, showing only а moderate deviation from theory. 


h. The variations in the magnitude of the moment derivative and 
in its phase angle, resulted in a trend toward instability at supercriti- 
cal Mach numbers. In particular, for airfoils of equal thickness the 
effect of an extreme forward location of maximum thickness was destabiliz- 
ing in that negative aerodynamic damping existed, implying the possibility 1 
of a single degree of freedom type of flutter. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., Mar. 24, 195} 
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TABLE I.- MEASURED FLUTTER DERIVATIVES FOR THE NACA 65A012 AIRFOIL 
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TABLE II.- MEASURED FLUTTER DERIVATIVES FOR THE NACA 65A008 AIRFOIL 
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TABLE III.- MEASURED FLUTTER DERIVATIVES FOR THE NACA 65A004 AIRFOIL 
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TABLE IV.- MEASURED FLUTTER DERIVATIVES FOR THE NACA 2-008 AIRFOIL 
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TABLE V.- MEASURED FLUTTER DERIVATIVES FOR THE NACA 877A008 AIRFOIL 
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Figure 1.- Section profiles and pressure-cell locations of models. 
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Figure 2.- View of test section with model in place and diagrammatic sketch of drive system. 
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Figure 3.- Typical oscillograph traces. 
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Figure 4- Results as a function of reduced frequency, 
k, for several Mach numbers for the reference 
model, МАСА 654008; @ „= O°. 
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Figure Š.- Typical results for reference model, NACA 
654008; 2,0% 
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Figure 6.- Variation of experimental results from theory for 
reference model, NACA 65A008, witha faired line to 


show the mean variation with Mach number; а= 0°. 
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Figure 7.- Concluded. 
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Figure 8.- Effect of airfoil thickness on lift derivatives. 
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Figure 8.- Concluded. 
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Figure 10- Comparison of steady and unsteady lift derivatives 
for airfoils with varying thickness. 
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Figure /O.- Concluded. 
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Figure /L- Results as a function of reduced frequency, 
k, for several Mach numbers for the reference 
model, МАСА 654008; @„= O°. 
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Figure /2 - Variation of experimental results from theory for 
reference model, NACA 654008, witha faired line to 
show the mean variation with Mach number; ат= 0°. 
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Figure /3.- Effect of airfoil thickness distribution on momen? 


derivatives. 
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Figure /4.- Effect of airfoil thickness on moment derivatives. 
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(а) Effect of airfoil thickness distribution. 
Figure /5.- Damping component of the moment derivatives. 
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(b) Effect of airfoil thickness, 
Figure /5.- Concluded. 
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